Confidence Limits for the Ratio of Two Binomial Proportions Based on Likelihood Scores: Non-Iterative Method

The interval estimation of the ratio of two binomial proportions based on the score statistic is superior over other methods. Iterative algorithms for calculating the approximate confidence interval have been provided by, e.g., KOOPMAN (1984, Biometrics 40:513–517) and GART and NAM (1988a, Biometrics 44:323–338). This note presents the analytical solutions for upper and lower confidence limits in a closed form and gives examples for numerical illustration. The non-iterative method is generally more desirable than the iterative method.     

Genetic structure and diversity of the Diachasmimorpha longicaudata species complex in Thailand: SSCP analysis of mitochondrial 16S rDNA and COI DNA sequences

Genetic variation among 20 populations of Diachasmimorpha longicaudata (Ashmead) (Hymenoptera: Braconidae) in Thailand was investigated using single strand conformation polymorphism (SSCP) analysis of mitochondrial DNA sequences. From a total 641 individual parasitoids, seven distinct haplotypes containing a total of 32 polymorphic sites were observed from cytochrome oxidase subunit I (COI) sequences along with five distinct haplotypes containing a total 16 polymorphic sites from 16S rDNA sequences. Values obtained through pairwise F_ST comparisons and analysis of molecular variance (AMOVA) indicated significant genetic differentiation among D. longicaudata populations in Thailand. Congruent relationships showing separation of these populations into three groups were obtained from Neighbor joining and Bayesian phylogenetic tree analyses along with the use of haplotype networks. This SSCP analysis of populations of the D. longicaudata species complex is the first report using molecular population genetic methods to analyze the structure of this parasitoid species in Thailand. This may provide useful information for release of parasitoid strains to maximize their benefit in biological control programs.     

Determination of protein and amino acid requirements of lactating sows using a population-based factorial approach

Determination of appropriate nutritional requirements is essential to optimize the productivity and longevity of lactating sows. The current recommendations for requirements do not consider the large variation between animals. Therefore, the aim of this study was to determine the amino acid recommendations for lactating sows using a stochastic modeling approach that integrates population variation and uncertainty of key parameters into establishing nutritional recommendations for lactating sows. The requirement for individual sows was calculated using a factorial approach by adding the requirement for maintenance and milk. The energy balance of the sows was either negative or zero depending on feed intake being a limiting factor. Some parameters in the model were sow-specific and others were population-specific, depending on state of knowledge. Each simulation was for 1000 sows repeated 100 times using Monte Carlo simulation techniques. BW, back fat thickness of the sow, litter size (LS), average litter gain (LG), dietary energy density and feed intake were inputs to the model. The model was tested using results from the literature, and the values were all within ±1 s.d. of the estimated requirements. Simulations were made for a group of low- (LS=10 (s.d.=1), LG=2 kg/day (s.d.=0.6)), medium- (LS=12 (s.d.=1), LG=2.5 kg/day (s.d.=0.6)) and high-producing (LS=14 (s.d.=1), LG=3.5 kg/day (s.d.=0.6)) sows, where the average requirement was the result. In another simulation, the requirements were estimated for each week of lactation. The results were given as the median and s.d. The average daily standardized ileal digestible (SID) protein and lysine requirements for low-, medium- and high-producing sows were 623 (CV=2.5%) and 45.1 (CV=4.8%); 765 (CV=4.9%) and 54.7 (CV=7.0%); and 996 (CV=8.5%) and 70.8 g/day (CV=9.6%), respectively. The SID protein and lysine requirements were lowest at week 1, intermediate at week 2 and 4 and the highest at week 3 of lactation. The model is a valuable tool to develop new feeding strategies by taking into account the variable requirement between groups of sows and changes during lactation. The inclusion of between-sow variation gives information on safety margins when developing new dietary recommendations of amino acids and protein for lactating sows.     

Partitioning of knee joint internal forces in gait is dictated by the knee adduction angle and not by the knee adduction moment

Medial knee osteoarthritis is a debilitating disease. Surgical and conservative interventions are performed to manage its progression via reduction of load on the medial compartment or equivalently its surrogate measure, the external adduction moment. However, some studies have questioned a correlation between the medial load and adduction moment. Using a musculoskeletal model of the lower extremity driven by kinematics–kinetics of asymptomatic subjects at gait midstance, we aim here to quantify the relative effects of changes in the knee adduction angle versus changes in the adduction moment on the joint response and medial/lateral load partitioning. The reference adduction rotation of 1.6° is altered by ±1.5° to 3.1° and 0.1° or the knee reference adduction moment of 17 N m is varied by ±50% to 25.5 N m and 8.5 N m. Quadriceps, hamstrings and tibiofemoral contact forces substantially increased as adduction angle dropped and diminished as it increased. The medial/lateral ratio of contact forces slightly altered by changes in the adduction moment but a larger adduction rotation hugely increased this ratio from 8.8 to a 90 while in contrast a smaller adduction rotation yielded a more uniform distribution. If the aim in an intervention is to diminish the medial contact force and medial/lateral load ratio, a drop of 1.5° in adduction angle is much more effective (causing respectively 12% and 80% decreases) than a reduction of 50% in the adduction moment (causing respectively 4% and 13% decreases). Substantial role of changes in adduction angle is due to the associated alterations in joint nonlinear passive resistance. These findings explain the poor correlation between knee adduction moment and tibiofemoral compartment loading during gait suggesting that the internal load partitioning is dictated by the joint adduction angle.     

Temporal changes of mechanical signals and extracellular composition in human intervertebral disc during degenerative progression

In this study, a three-dimensional finite element model was used to investigate the changes in tissue composition and mechanical signals within human lumbar intervertebral disc during the degenerative progression. This model was developed based on the cell-activity coupled mechano-electrochemical mixture theory. The disc degeneration was simulated by lowering nutrition levels at disc boundaries, and the temporal and spatial distributions of the fixed charge density, water content, fluid pressure, Von Mises stress, and disc deformation were analyzed. Results showed that fixed charge density, fluid pressure, and water content decreased significantly in the nucleus pulposus (NP) and the inner to middle annulus fibrosus (AF) regions of the degenerative disc. It was found that, with degenerative progression, the Von Mises stress (relative to that at healthy state) increased within the disc, with a larger increase in the outer AF region. Both the disc volume and height decreased with the degenerative progression. The predicted results of fluid pressure change in the NP were consistent with experimental findings in the literature. The knowledge of the variations of temporal and spatial distributions of composition and mechanical signals within the human IVDs provide a better understanding of the progression of disc degeneration.     

On partial local smoothing rules for curve estimation

We compare the performances of local and global rules for smoothingparameter choice, in terms of asymptotic mean squared errorsof the resulting estimators. In some instances there is surprisinglylittle to choose between local and global approaches; our analysisidentifies contexts where the differences are small or large.This work motivates development of smoothing rules that forma ‘half-way house’ between local and global smoothing.There, interpolation provides a basis for partial local smoothing.A key result shows that interpolation on even a coarse gridcan produce a very good approximation to full local smoothing.Our theoretical and numerical results lead us to suggest linearinterpolation of a bandwidth obtained by integral approximationson discrete intervals.